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Abstract 

We study a CP-violating triple product correlation that occurs in the decay of a neutral Higgs 
boson into tt pairs when the Higgs boson does not have a definite CP nature. We consider the 
H ^ tt decay channel as well as the gluon fusion process gg ^ H ^ ti. The asymmetry in Higgs 
decay, normalized to the H ^ tt width, can reach the 6% percent level. In the gluon fusion process 
the corresponding normalized asymmetry is smaller by an order of magnitude. We present a crude 
estimate of this observable at the LHC. 
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I. INTRODUCTION 


In models with more than one Higgs doublet it is possible to have CP violation in the 
coupling of neutral Higgs bosons to fermions. Such CP violation manifests itself in the form 
of a physical neutral Higgs that has both scalar and pseudo-scalar couplings. The couplings 
to the top-quark would then be given by the Lagrangian 

C = ——Ht{A + iB'y^)t, (1) 

V 

where A,B are real and H = 1, H = 0 corresponds to the standard model with one Higgs 
doublet. Exploring the possible existence of these couplings will be an important task for 
the LHC. 

The determination of the CP nature of a neutral Higgs boson at the LHC has been 
studied before in the literature. Gunion and He considered the production of a light 
Higgs boson in association with a ti pair Ending that the spin-averaged cross-section does 
not contain terms proportional to AB. These terms are also absent in the spin-averaged 
cross-section for the process gg —>■ ti. However, the weak decay of the top quark makes it 
possible to study spin correlations because it analyzes the top polarization P. The spin 
correlations manifest themselves in the angular distributions of the decay products of the t 
and i. These top-quark spin correlations have been the subject of much attention [^, and 
in our case they provide an avenue to explore the product of couplings AB. 

A particular example of a CP violating observable that is proportional to AB is the 
T-odd triple product correlation 0,1^ that arises for the decay chain H —>■ ti ^ bbW~^W~: 

^{PuPuPb.Pb) = ^{Pw+,Pw-,Pb,Pb) = VsPb ■ {pt X Pb) (2) 

where the last equality follows in the ti (parton) center of mass frame. In this frame, the 
difierent momenta transform under CP as, 

Pt Oil -Pi = Pt , Pb ^ -Pb- (3) 

From these transformation properties it follows that this triple product correlation is a 
CP-odd observable: 


Pt ■ {Pb X Pb) QE, -pt ■ {Pb X Pb)- (4) 

The correlation can then be used to signal and measure CP violation in Higgs decay as well 
as in the gluon fusion process gg ^ ti ^ bbW^W~. At the LHC the pp initial state is not a 
CP eigenstate and the signal could receive contributions from CP conserving backgrounds. 
This has been discussed previously for diEerent CP-odd observables considered in the lit¬ 
erature j^. In particular, the CP odd correlation involving the ti spin and momenta that 
underlies our asymmetry has been previously discussed by Bernreuther and Brandenburg 

i- 

In this paper we calculate the CP violating observable Eq.|21for both the Higgs decay and 
gluon fusion processes. We also present a crude estimate for the LHC but leave the details 
of a more realistic simulation, including the study of background, to a future publication. 
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II. HIGGS DECAY 


We begin by considering the process H ^ tt ^ bbW~^W for the interaction given by 
Eq. HJ The total Higgs decay width into a tt pair is then given by 


T{H ^tt) = N, 
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The decay distribution for this process contains a triple product correlation of the form 
Eq. 121 This CP-violating correlation can be obtained in a straightforward manner; after 
summing over the spin of the b and b as well as the polarizations (collectively denoted 
by s) and the quark color, we hnd: 





( 6 ) 


For the intermediate tt states we have used the narrow width approximation. We assume 
that the f-decay occurs as in the standard model and that = 1, therefore 
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To proceed we write the four body phase space choosing as independent variables the t 
angles in the Higgs rest frame, the b angles in the t rest frame and the b angles in the i rest 
frame. The integration range for this choice of angles is not constrained and we write 


= 4m, - 4 ;;^^^*^ 


( 8 ) 


where * denotes expressions evaluated in the rest frame of t and ** expressions evaluated in 
the rest frame of t. 

In the Higgs rest frame the correlation can be simplihed further as illustrated in Figure ^ 
We use the t-quark momentum in this frame to dehne the z-axis and the t bW^ decay 
plane to dehne the x — z plane. The i bW~ decay plane then forms an angle (pf* with 
respect to the f-decay plane. With these choices we hnd 


4Pt,Pi,Pb,Pb) = MhPi ■ {pb X pt) 

= -MulPtlPbxPby (9) 

We then note that pbx and piy are invariant under a boost to the t and i rest-frames respec¬ 
tively, and use this to calculate them in those frames. This leaves us with a correlation of 
the form 


^{PuPuPb.Pb) = M, 



sin 6t sin OP 


sm 0T. 


h ■ 


( 10 ) 
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FIG. 1: Kinematics for the reaction H ^ tt ^ bbW~^W~. The t-quark momentum, pt, defines the 
z-axis. The t-quark decay defines the x — z plane and the t decay plane forms an angle (j)^* with 
respect to the t decay plane. 


After this, the integration over phase space is straight-forward. Dehning an integrated CP 
odd observable T: 


f = J dT{H ^ tt ^ bbW^W ) sign{e{pt,PuPb,Pb)) 
= J dT{H bbW+W-) sign{sm(j)l*) 


we arrive at our normalized CP asymmetry, 
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This normalized observable corresponds to the simple counting asymmetry 


^eventsiPb ' i.Pb ^ Pt) ^0) ^events i^Pb ' i.Pb ^ Pt) ^ 0) 
^eventsiPb ' i.Pb ^ Pt) ^ 0) -|- Neventsi^Pb ' {pb ^ Pt) ^ 0) 


( 11 ) 


( 12 ) 


(13) 


in the Higgs rest frame. 

Weinberg has shown that there are unitarity constraints on the size of the CP violating 
couplings of the type of Eq. [TJ For example, in the models discussed by Weinberg , the 
product AB from Eq. ^corresponds to what he calls Im Z 2 . Assuming that the lightest 
neutral Higgs eigenstate dominates and that the different scalar vacuum expectation values 
are comparable, Eqs. 49 or 52 of Ref. then imply. 




( 14 ) 
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Ref. also shows how to construct models that reach this bound For our numerical 
estimates of the CP-odd signal, Acp, we will choose A = B = 1/a/2 in accordance with 

Eq. [m 

Numerically this result is shown in Figure El The asymmetry can have either sign de- 



Mjj(GeV) 


FIG. 2: Normalized CP violating asymmetry in Higgs decay as given by Eo. 1121 We use A = B = 
1 /\/2 from Eq. El 

pending on the relative sign of A and B. It could also be quite large at over 6%, and is 
nearly independent of the Higgs mass (as long as the Higgs is sufficiently heavy to decay 
into ti). 

^ For comparison, the first paper in Ref. Q, also uses Eq. d and the papers in Ref. 0 use the equivalent 
of HR ~ 1 
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III. CP-ODD ASYMMETRY IN GLUON FUSION 


In hadron colliders the Higgs boson is produced by gluon fusion through an intermediate 
top-quark loop. This leads us to consider the combined Higgs production and decay through 
the reaction gg ti. For this source of ti pairs there are three relevant diagrams as shown 
in Figure El There is the resonant ti production shown in the first diagram, and there are 
additional t and u channel diagrams in which the gluons couple directly to the ti pairs which 
can interfere with the Higgs exchange diagram. The complete CP violating asymmetry in 



FIG. 3: Three diagrams responsible for CP asymmetry in top-quark pair production 

this case receives two contributions. The hrst contribution arises from the square of the 
hrst diagram and corresponds to the asymmetry in Higgs decay modihed by factors for the 
Higgs production vertex. The second contribution arises from the interference between the 
Higgs exchange diagram and the t and u channel diagrams. Note that there is an additional 
s-channel diagram in which a gluon couples to the ti pair. However, this diagram does not 
interfere with the Higgs exchange diagram and does not contribute to the CP asymmetry. 

The effective coupling between the Higgs boson and gluons at one-loop can be written in 
the form 


= - [Pa(s)(-gi ■ q29^,^ + q 2 f,qiu) + Fb(s)e(/r, u, qi, ga)] • (15) 


The two form factors Fa{s) and Fb{s) have been obtained previously in the literature. For 
the kinematic regime in which the Higgs boson is heavier than a ti pair they are given by 
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This is the case where ^/s is large enough to produce top-quark pairs, and therefore the form 
factors are complex. CP odd but naive-T even observables (some of which are discussed in 
Ref. j^) require a “strong phase” and are thus proportional to the absorptive part of these 
form factors. Our UP-odd triple product correlation, on the other hand, does not require a 
“strong phase” and will be proportional to the real part of these form factors 0]. 
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A convenient way to calculate the CP asymmetry is to consider the process as in Figure 0] 
in the parton CM frame and use helicity amplitudes for the most part. The top-quark pair 
production by the three diagrams in Figure 121 is represented by Fp in Figured The t and t 
decays into hW are represented by Fp p. It follows that the amplitude can be written as 



FIG. 4: Decomposition of tt production and decay vertices with helicity amplitudes. 


Ub^pHit + mt)Tp{-i)t + m^)F^n5e*^(fF+)e*^(fF ) 
(P? - rn^t){Pi - mj) 


(18) 


To construct helicity amplitudes we now replace the numerator of the top-quark (and anti¬ 
top-quark) propagator with a sum over polarizations. We work within the narrow width 
approximation for the t and i decays and, therefore, these polarization sums refer to on-shell 
ti states. As 
in which the 
replacement: 


We now sum 

■ hu'7MA7v(l - 75)MiA (^- 9 ^’' + ) 

• D-^7m'I^67i^'( 1 - 75)wtV (^- 9 ^'"'' + j 

■ UtxTpVia hiV7°rJ,7°Mu' (20) 


we will see, the CP odd observable arises from the interference of amplitudes 
intermediate states have different helicities. Our amplitude becomes after this 


M = 


Ub'^DiT.xUtxUt\)Tp{J2a'ViaVia)'^f)Vie*^{W+)e*''{W ) 


{Pt - m^){p‘l 


ml) 


(19) 


over the and the W polarization as well as the b and the b spin to obtain: 
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and we have already used the standard model decay vertex 

The last line in Eq. 1201 corresponding to the production vertex, is more easily computed 
in the helicity basis and we write it as 

Ut\T T X! -^2, A, (T) A^p*(Ai, A 2 , A', Cx') (21) 

Ai,A2 

where Mp represents the helicity amplitudes for gg —>■ ti and the factor 1/4 arises when we 
average over the gluon helicities. These helicity amplitudes are summarized in Tabled 
The second line of Eq. 1201 corresponding to the f-decay, can be simplihed as follows 

- 75)utx 

^ ““ (“‘ (ii; ■ 0 ^ ■ i|:j) 

The f-decay factor simplihes in a similar manner 

^ (™‘ * (" ■ i|:)) 

In these last two factors, Eqs. [2211201 only the terms proportional to ji;, and contribute 
to the CP-odd asymmetry (as otherwise there are not enough independent four vectors to 
form the triple product correlation). The surviving terms are also calculated as a function 
of the t and i helicities and we use the dehnitions. 


7)(A',A) = Utx'M^ - 75)utx 

Ti{a,a') = hto|^b(l - 75)^'^v• (24) 


These factors are given explicitly in Table iTTl With all this we arrive at the result. 
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MriXi, A2, a, a)Mp*(Ai, A 2 , A', X)Tt{a, a') 


(25) 


To evaluate this expression we use Table d and Table O in the appendix. They are 
calculated in the parton (gg) center of mass frame where the t and b four-momenta are 
parametrized by 


/ s 

Pt = “^(1) /5sin 6*, 0, /3 cos 6*), pb = Pfe(l, sin6*b cos 0b, sin sin 0b, cos^b), (26) 

and similarly for the h. (3 = — druf/ s is the usual kinematic factor. We emphasize that 

these angles correspond to the parton CM frame, unlike the angles used in Section II, which 
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were given in the t and t rest frames. In terms of these kinematic variables the desired triple 
product correlation reads 

e{pw+,Pw-,Pb,Pb) = ^{PuPhPb.Pb) = ~2^bEi 

(sin( 6 ') cos( 6 *b) sin( 6 ';,) sin( 0 b) — cos( 6 ') cos( 05 ) sin( 6 'b) sin( 6 ' 5 ) sin(0fe) 

— sin( 6 ') cos( 6 *b) sin( 6 *b) sin(05) + cos( 6 *) cos(0fe) sin( 6 *f,) sin( 6 * 5 ) sin(0^)). (27) 



The hrst contribution to this correlation, from the square of the Higgs exchange diagram, 
can be extracted from Eq. |2S1 Its corresponding color factor is 5a,b^a,bNcI9>^ = iV^/S and 
we End: 
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(28) 


We have written Ea. 1281 explicitly in terms of Eq. IHl This serves as a check of the calculation 
since this result can also be obtained from multiplying Eq. El with appropriate factors for the 
Higgs propagator and the Higgs production vertex. The resulting correlation is not linear 
in any absorptive phase as expected. The subscript tpc indicates that we have only retained 
those terms that give rise to the CP violating correlation and dropped all others. 

The second contribution to the asymmetry results from the interference of either the 
t-channel or w-channel gg —^ tt diagram with the s-channel Higgs diagram. It can also be 
extracted from Eq. |2Slbut this time the color factor is equal to 6a,b'^t^(T^T^)/8‘^ = 1/16. 
The desired correlation is found to be proportional to the coupling constant combination 
63fJ[(Fa)] + a3fJ[(Fft)], once again no absorptive phase is required to generate the CP violating 
signal. Our result for this contribution to the asymmetry in the parton center of mass frame 
is 
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€{Pt,Pi,Pb,Pb)SiPt - mt)^iPl - ^t) 
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Tt [1 + (2 - Tt)^f{Tt)]] 

(s — Mfj) (^sin^ 0 + Tt cos^ O'j 

(29) 


The complete asymmetry in gluon fusion is the sum of Eq. |2Hl and Eq. (221 In the tt center 
of mass frame it corresponds, once again, to the simple counting asymmetry of Eq. ^1 Of 
course, these results need to be folded with appropriate parton distribution functions and 
evaluated in the lab frame to arrive at physical observables for hadron colliders. We do this 
for the case of the LHC in the next section. 


IV. RESULTS FOR THE LHC 

At the LHC there are three main mechanisms for tt production: gluon fusion, resonant 
production via an intermediate Higgs boson and light qq annihilation. We reproduce in 
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Figure El the relative magnitudes of these mechanisms with Higgs couplings given by Eq. ^ 
with A = B = I/V 2 per Eq. CH as a function of the tt pair invariant mass The resonant 
production is shown for two values of the Higgs mass. This hgure illustrates the well known 



M„(GeV) 


FIG. 5: ti production mechanisms at the LHC. The larger cross-section indicated by the dashed 
line corresponds to gluon fusion. The dotted line indicates light qq annihilation. We also show the 
resonant Higgs production for two different values of the Higgs mass. 

result that the gluon fusion contribution dominates for most of the rriu range (except very 
^ Throughout this section we will use the CTEQ6M parton distribution functions and y/S = 14 TeV. 
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close to the Higgs resonance where the two mechanisms are comparable) . The light qq 
production mechanism is about an order of magnitude smaller for Mu ^ 750 GeV. 

We hrst study the CP violating asymmetry generated by the resonant production of a 
Higgs at the LHC. This contribution to the total asymmetry is obtained from Eq. |2H1 after 
we use the narrow width approximation for the Higgs propagator. In the lab frame the 
correlation cannot be rewritten as the simpler triple product pi ■ {pb x p)) of Eq. |21 Instead, 
the physical asymmetry is now dehned by 


Nevents {^{Pu Ph Pb, Pb) > 0) " Neyents {^{Pu Ph Pb, Pb) < 0) 

^events : Pt5 Pb) ^ 0 ) “f* -/Vguenbs 7 Pb? Pb) 0 ) 


(30) 


and requires full reconstruction of four four-momenta. 

It is possible to construct a different GP-odd correlation that generalizes the triple prod¬ 
uct of Eq. 121 to the lab frame and that only requires the reconstruction of the directions of 
the relevant momenta. One such example is 


^events {.{Pt Pt^ ' {,Pb ^ Pb) ^ 0) ^events {.{Pt Pt) ' {Pb ^ Pb) ^ 0) 
Nevents {{Pt “ Pt) ' (Pb X Pg) > 0) + Nevents {{Pi “ Pt) ' {Pb X Pb) < 0) ‘ 


(31) 


In Figure El we show both Acp (solid line) and Acp (dotted line) as a function of the 
Higgs mass. We see that the production factors in Eq. 12H1 cancel out and Acp is the same 
as calculated for Higgs decay (in the narrow-width s-channel approximation) and shown in 
Figure 121 On the other hand, there is a dilution in the asymmetry of about 30% when using 
the form Acp- 

The complete asymmetry at the LHC is obtained by adding Eqsl2Hl and 1221 (and folding 
in the corresponding parton distribution functions) The normalized asymmetry is shown 
in FigureEl We present results for Acp (solid line) and Acp (dotted line). Both of them are 
signihcantly smaller than in Higgs decay because the total cross-section for tt production 
(which appears in the denominator of Eq. EOl and of Eq. I!?T1) is much larger than the Higgs 
production cross-section. The asymmetry does not increase as much as the cross-section 
because the gluon fusion mechanism only contributes to it through interference with the 
Higgs resonance. We hnd once more that Acp is smaller than Acp by about 30%. Notice 
that the light qq tt process can only contribute a very small amount to the asymmetry. 
This is because it occurs dominantly through one gluon exchange which does not interfere 
with the resonant Higgs production channel. 


V. SUMMARY AND CONCLUSIONS 

We have presented a new CP odd observable in the form of a triple product correlation 
that is sensitive to a mixture of scalar and pseudoscalar couplings of a Higgs boson to the 
top-quark. We have evaluated this observable in simple Higgs decay as well as in gluon 


^ There is also a small contribution to the asymmetry from light qq Higgs production which we neglect in 
our discussion. 
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FIG. 6: CP Asymmetry for s-channel Higgs production and subsequent decay at the LHC. The 
solid and dotted curves correspond to the asymmetries Acp, Eq. EH and Acp, Eq. Irrespectively. 

fusion production of ti pairs and illustrated that the resulting asymmetry can occur at the 
6% percent level if the intrinsic CP violation is near its unitarity bound. To reconstruct the 
asymmetry it is necessary to observe the t (or t) decay plane as well as the direction of the 
b (or b) momentum. 

The observable originates in a tt spin correlation that is analyzed by the weak decay of 
the top-quark. For this reason it will only apply to very heavy Higgs bosons that can decay 
into ti pairs. For lighter Higgs bosons it is possible to construct a similar observable in the 
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FIG. 7: Normalized CP violating asymmetry at the LHC for A,B = The solid and dotted 

curves correspond to the asymmetries Acp, Eq. EOl and Acp, Eq. respectively. In this case 
both contributions to the asymmetry are included and they are normalized to the total pp —> tt 
cross-section. 

H —>■ T^T~ channel. Althongh this extension is straightforward, it is not clear that it will 
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be possible to identify pairs and their decay planes at high energy colliders. 

We have estimated the size of the asymmetry in models with multiple Higgs bosons as the 
source of CP violation. However, the asymmetries Aqp from Eq. ISOl and Acp from Eq. 
are more general than this and can be generated by other mechanisms of CP violation. In 
this sense they constitute a new tool to search for CP violation at future colliders. 

We have illustrated a preliminary application of our result to pp collisions at the LHC. In 
this case reconstruction of the signal Acp in the lab frame is harder and there is a dilution 
if one uses the simpler signal Acp instead. For leptonic W decays the direction of the 
t{t) cannot be reconstructed. It is then necessary to dehne a different asymmetry. One 
possibility is to replace the t{t) momenta with the momenta of the leptons from the W 
decay. Denoting by an or a pA pair one could use, for example. 


^events ((Pr+ Pi~^ ' (Pfe ^ Vb) ^ 0) ^events ((Pr+ ) ' (Pfe ^ Pb) ^ 0) 

^events ((P£+ P£“) ' (Pb ^ Pb) ^ 0) p N events ((P£+ P£~) ' {.Pb ^ Pb) ^ 0) 


(32) 


However, we expect a further dilution of the signal in this case. A further complication at 
the LHC is that the initial pp state does not transform into itself under CP conjugation 
and therefore, it is possible to have CP conserving backgrounds to these asymmetries. A 
discussion of these backgrounds, as well as a more detailed simulation of the asymmetry 
are necessary for a realistic evaluation of the usefulness of this new observable at the LHC. 
We leave these considerations to a future publication. There exist several papers in the 
literature that develop CP asymmetries in similar decay chains to the one considered here in 
the context of or 77 colliders [n|. The papers of Ref. 0 * hrst discussed the underlying 
correlation involving the spin and momenta of the tt pair. They developed the asymmetry 
in a slightly different manner that involves the momenta of the secondary leptons. It is 
difficult to compare their results with our form for the asymmetry without a more detailed 
numerical simulation. 
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TABLE 1: t + u channel gg —> tt and gg ^ H ^ tt helicity amplitudes. The former should 
be multiplied by g^ and the latter by {gmt/2Mw)- Color factors have not been included and 
the kinematic factors refer to the t-quark in the parton CM frame: Tt = ^rn^js^ f3 = -y/1 — Tt, 
7 = (1-/32)-V2. 
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APPENDIX A: TABLES OF HELICITY FACTORS 

In this appendix we present our helicity factors. They are calculated in the parton center 
of mass frame with the first parton (g) defining the z-axis and 6 being the polar angle of 
the t-quark momentum. These expressions are used in this paper for pp collisions at the 
LHC, therefore terms that are odd under the interchange z —z will not contribute to the 
differential cross-section after the parton distribution functions are folded in. 

First we present the helicity amplitudes for gg —> ti. The kinematic factors in Table Q] 
refer to the top-quark. Only the first four rows of the table contribute to the asymmetry 
which requires interference with the s-channel Higgs exchange amplitude. The remaining 
rows are quoted for completeness. 

Finally we present in Table ITTl the t and t decay factors as defined in Eqs. In addition 
to the top-quark kinematic factors, the results contain the parton center of mass polar and 
azimuthal angles for the b and b momenta as defined in Eq. |2S1 Notice that the factors that 
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TABLE II: t and t decay factors as defined in the text. The first column should be multiplied by 
Ef, and the second column by Eg, the b and b energies in the parton center of mass frame. 


A, A' 


Tt{\,\') 


■ 05(1 -|- / 3)(1 — COS 0 COS 6*6 — sin6* sin6*6 cos pb) 

■0s(l — / 3)(1 — cos 6* cos 6*5 — sin 0 sin 6*5 cos pi) 

+ 

2mt{— sin 6* cos 9 b -f cos 6* sin 6b cos pb — i sin 6b sin pb) 

2mt{— sin 6* cos 61 + cos 6* sin 6^ cos Pb + i sin 6^ sin pp) 

+ ) “ 

2mt{— sin 6* cos 6b -f cos 6* sin 6b cos pb + i sin 6b sin pb) 

2mt{— sin 6* cos 6^ -|- cos 6* sin 6^ cos p^ — i sin 6^ sin pp) 

+ ) + 

■05(1 — / 3)(1 -|- cos 6* cos 6*6 -|- sin 6* sin 6*6 cos pb) 

■ 0 s(l -|- (5){1 + cos 6* cos 6*5 -|- sin 0 sin6*5 cos pp) 


generate the triple product correlation, sin 6b sin 0 ^ and sin 6*5 sin (pi only occur in the factors 
with “mixed” helicities. That is, the correlation arises from the interference of different 
helicity amplitudes of the intermediate top-quark. 
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